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Abstract. The study of the gamma-ray radiation produced by cosmic rays that escape their 

accelerators is of paramount importance for (at least) two reasons: first, the detection of those 

gamma-ray photons can serve to identify the sources of cosmic rays and, second, the charac- 

£/■) teristics of that radiation give us constraints on the way in which cosmic rays propagate in the 

r- h interstellar medium. This paper reviews the present status of the field. 
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1. Introduction 



(f~) The galactic disk is the sky's most prominent source of gamma rays in the GeV energy 

domain. Discovered in the late sixties (Clark et al. 1968), the galactic gamma-ray diffuse 

emission was soon interpreted as the decay of neutral pions (tt°) produced by cosmic 

rays (CRs) interacting with the interstellar gas (Steckcr 1969). This confirmed a scenario 

first proposed by Hayakawa in 1952 (for early reviews see Fazio 1967 and Ginzburg & 

rS^ Syrovatskii 1964). During the past decades the diffuse gamma-ray emission from the 

Wh galactic disk has been observed with constantly increasing accuracy by several space 

instruments (Fichtel et al. 1975, Mayer-Hasselwander et al. 1982, Hunter et al. 1997, 

Ackermann et al. 2012a), and the spatial distribution of CRs in the Galaxy could be 

extracted from such observations (Bloemen 1989 and references therein, Bertsch et al. 

i i 1993, Strong et al. 2004, Ackermann et al. 2012). It turned out that the distribution of 

CRs in the Galaxy is quite uniform and that, as an order of magnitude, the intensity of 

CRs measured in the solar system is representative of the intensity anywhere else in the 

galactic disk. This roughly uniform background is often referred to as the cosmic ray sea. 

j^j The 7r°-decay gamma-ray emissivity q~(> 100 MeV) of the local atomic gas has been 

^H measured by the Fermi LAT and is equal to <j°(> 100 MeV)/4-7r = 1.6x 10~ 26 ph/s/sr/H- 

i-H atom (Abdo et al. 2009). Thus, an estimate of the intensity of the galactic diffuse emission 

_il from a specific direction in the sky can be obtained by integrating the gamma-ray emis- 

Q sivity of the gas q 7 (> 100 MeV, I) w q°(> 100 MeV) along the line of sight I. This 

CO straightforward procedure gives (e.g. Aharonian 2004): 

(? 7 (>100MeV,Z) _, Cv ,, n -4/ N H \ ,,, 2 



> J 7 (> 100 MeV) = / dZ ^ ^ n gas w 1.6xl0" 4 1A?2 - „ ph/s/cmVsr 

.^_i J 47T \10^ cm l ) 

X . , (°) 

where Nh is a typical gas column density in the galactic disk. 

In 1973 Black & Fazio suggested that the gamma-ray emission from individual massive 
molecular clouds (MCs) might be visible above the diffuse galactic emission. The gamma 
ray flux from a MC of mass M c i, distance d and embedded in the CR sea is: 

<7°(> 100 MeV) (M cl /ra v ) „ / M* \ 

F 7 (> 100 MeV) w ^ ; V ' pJ w 2 x 10" 7 -^- ph/cm 2 /s (1.2) 
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where M5 is the mass of the cloud in units of 1O 5 M , dk pc its distance in kiloparsecs, 
and m p the proton mass. The detectability of massive MCs above the diffuse galactic 



emission follows from Equations (1.1) and (1.2 1 and from the fact that, for a typical 
cloud density of n c i w 1000 n 3 cm -3 , the radius of the cloud is R c i « 10 (M 5 /n 3 ) 1 / 3 pc 
and its angular extension is fl c i « 10~ 4 (M 5 /7i3) 2 / 3 dr ' sr. 



Black & Fazio (1973) proposed that the masses of the MCs of known distance could be 
derived from the strength of their gamma-ray flux , under the assumption that the inten- 
sity of CRs is known throughout the Galaxy (i.e. the assumption of a uniform CR sea). 
Such an approach was used extensively to calibrate the methods for the determination of 
the mass of molecular and atomic gas in the Galaxy and also to check the (rough) spatial 
homogeneity of the intensity of CRs in the Galaxy over large spatial scales (e.g. Car- 
aveo et al. 1980, Bloemen et al. 1984, Hunter et al. 1994, Digel et al. 1996, 1999, Abdo 
et al. 2010a, Ackermann et al. 2011, 2012b, 2012c, 2012d). However, the assumption of 
a uniform sea of CRs that permeates the whole Galaxy might be inappropriate in some 
circumstances, especially on small spatial scales (e.g. in the vicinity of CR sources). In 
these circumstances, the reasoning of Black & Fazio can be reversed and, if an estimate 
of the mass of the cloud is available (see e.g. Hartquist 1983 for caveats), one can use 
the gamma-ray observations of MCs to probe variations of the intensity of CRs in the 
Galaxy (Issa & Wolfendale 1981, Morfill et al. 1981, Aharonian 1991, 2001, Casanova 
et al. 2010). In this context, MCs serve as cosmic ray barometers, because from their 
gamma-ray flux it is possible to infer the intensity (and thus pressure) of the CRs. 

The fact that the study of MCs could help in solving the problem of the origin of 
CRs became evident when it was realized that an association between MCs and CR 
sources is indeed to be expected (Montmerle 1979). This is because CRs are believed to 
be accelerated at supernova remnant (SNR) shocks (e.g. Hillas 2005) and thus are likely 
to be produced in star forming environments (where core-collapse supernovae explode), 



which are in turn expected to host massive MCs. In a seminal paper, Montmerle (1979) 
described a scenario in which CRs, after being accelerated at SNRs, could escape the ac- 
celeration site and diffuse to a nearby MC and produce there gamma rays via interactions 
with the gas. Due to the presence of these runaway CRs, the CR intensity is expected to 
be strongly enhanced in the vicinity of SNRs. In the same way, the gamma-ray emission 
from a MC illuminated by the runaway CRs will be much larger than the one derived 



in Equation (1.2 1, which refers to a MC embedded in the CR sea. It follows that the 
detection of gamma rays from MCs located in the vicinity of SNRs might constitute an 
hint for the fact that the nearby SNR is (or was, in the past) acting as a CR accelerator. 

In this context, gamma-ray observations performed in the TeV domain are of great 
relevance, because: i) the expected gamma-ray emission for a MC illuminated by runaway 
CRs is expected to have a TeV flux which is within the reach of current Cherenkov 
instruments (Aharonian 1991, Aharonian & Atoyan 1996, Gabici et al. 2009), and ii) in 
the TeV energy domain the contribution from the CR sea (which has a steep spectrum) 
to the gamma-ray emission is virtually negligible, and thus any detection of MCs has 
necessarily to be interpreted as an excess of CRs above the sea at the location of the MC 



(Aharonian 1991). A spectacular example of this fact can be found in Aharonian et al. 
(2006), where the detection of a diffuse emission of TeV gamma rays was reported from 



a very massive MC complex located in the galactic centre region. These observations 
revealed an excess above the CR sea of a factor of w 4... 10 at TeV energies, and also 
an harder spectrum of CRs there (with slope ss 2.3). The excess indicates that a source 
(or more sources) of CRs might be present in the region (remarkably, runaway CRs from 
only one SNR would suffice to explain the observed gamma-ray emission). 
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A phenomenological description of the propagation of the CRs after their escape from 
the acceleration site has been developed in a pioneering paper by |Aharonian fc Atoyan| 
(1996)| who also discussed the expected radiative signatures (especially in gamma rays) 



due to the interactions of CRs in the ambient gas. Aharonian & Atoyan considered 
an isotropic and spatially homogeneous diffusion coefficient for CRs and computed the 
expected gamma-ray emission from the MC, and stressed the fact that the properties of 
such emission (intensity, spectral shape, duration in time, etc.) strongly depend on the 
value of the diffusion coefficient. This opens the possibility to constrain, from gamma-ray 
observations, the diffusion coefficient of CRs in the vicinity of their sources. This fact 
has a tremendous importance for CR studies, given that the diffusion coefficient is a very 
poorly determined quantity (both from an observational and theoretical point of view) . 
The multiwavclcngth emission resulting from the interactions of runaway CRs in a MC 
has been computed, for the specific case in which a SNR accelerates the CRs, by |GabicT 



& Aharonian (2007) and|Gabici et al. (2009) These studies have then been applied to 



specific situations in order to obtain constraints on the particle diffusion coefficient in 
the vicinity of SNRs (Gabici et al. 2010, Nava & Gabici 2013). Moreover, the TeV diffuse 
emission resulting from the interactions of runaway CRs in the diffuse interstellar medium 
(i.e. in the absence of a massive MC) has been predicted and found to be within the reach 
of future ground based instruments such as the Cherenkov Telescope Array, which is thus 
expected to play a crucial role in proving (or falsifying) the SNR paradigm for the origin 
of galactic CRs (Casanova et al. 2010, Acero et al. 2013). This paper is intended as a 



short review of these results. A more extended discussion can be found in Gabici (2013) 
Though the effectiveness of the penetration of CRs into MCs still remains an open 
issue (Skilling & Strong 1976, Cesarsky & Volk 1978, Morfill 1982, Zweibel & Shull 1982, 
Everett & Zweibel 2011) I will assume in the following a full, unimpeded penetration of 
CRs in MCs. There is little doubt that TeV CRs can penetrate MCs (Gabici et al. 2007) 
and the penetration of GeV particles seems to be supported by gamma-ray observations, 
from the early ones by Lebrun fc Paul (1978)|to the ones by Ackermann et al. (2012b) For 



lower energies this remains an open issue, but this should not affect our considerations. 
Here, I will not discuss the (yet not clear) way in which CRs escape their sources 
(see Gabici 2011 and references therein) nor the case in which the gamma-ray emission 
from the MC is the result of the interaction between the cloud itself and the SNR shock 
(numerous papers can be found in the literature, including: Blandford & Cowie 1982, 
Aharonian et al. 1996, Gaisser et al. 1998, Bykov et al. 2000, Fatuzzo & Melia 2005, 
Uchiyama et al. 2010, Malkov et al. 2011, Inoue et al. 2012, Fang & Zhang 2013) 

2. Escape of cosmic rays from supernova remnants: isotropic diffusion 

As an illustrative example I consider here the case of the SNR W28, which is an aged 
remnant (t age w 4 x 10 4 yr) located in the vicinity of three massive MCs (of total mass 
« 1O 5 M ). Gamma-ray emission has been detected from the MCs in both the GeV and 
TeV domain (Aharonian et al. 2008, Abdo et al. 2010b, Giuliani et al. 2010). Since most 
of the gamma-ray emission clearly comes from outside of the SNR shell, it seems natural 
to interpret it as the result of the interactions in the MCs of CRs that escaped the SNR 
(Fujita et al. 2009, Gabici et al. 2010, Li & Chen 2010, Ohira et al. 2011, Yan et al. 
2012). This supports the idea that W28 was, in the past, an accelerator of CRs. 

I provide now a simple argument to show how one can attempt to constrain the diffusion 
coefficient in the vicinity of the SNR W28 by using the above mentioned gamma-ray 
observations, especially the ones performed by H.E.S.S. in the TeV domain. The time 
elapsed since CRs with a given energy E escaped the SNR can be written as: tdiff = 
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tage — tesc, where t esc (E) is the age of the SNR when CRs of energy E were released. This 
is a time dependent quantity, since the highest energy CRs are believed to be released 
first, and CRs with lower and lower energy are gradually released at later times (Gabici 
2011 and references therein). However, for CRs with energies above 1 TeV (the ones 
responsible for the very high energy gamma-ray emission) one can assume t esc << t age 
(i.e. high energy CRs are released when the SNR is much younger than it is now) and 
thus tdiff ~ t age . At time t age CRs have diffused over a distance Rd ~ a/4 D t age . 
Within the diffusion radius Rd the spatial distribution of CRs, fcR, is roughly constant, 
and proportional to ijEsn / '-R^, where E$n is the supernova explosion energy and T] is 
the fraction of such energy converted into CRs. On the other hand, the observed gamma 
ray flux from each one of the three MCs detected in gamma rays is: F 1 ex fcRM c i/d 2 , 
where M c i is the mass of the MC and d is the distance of the system. Note that in 
this expression F 1 is calculated at a photon energy E 7 , while Jcr is calculated at a CR 
energy Ecr ~ 10 x E 7 , to account for the inelasticity of proton-proton interactions. By 
using the definitions of Jcr an d Rd one can finally write the approximate equation, valid 
within a distance Rd from the SNR: 

p V E SN / M cl 

7 K (x D gal t age y/ 2 \ d 2 

Estimates can be obtained for all the physical quantities in the equation except for the 
CR acceleration efficiency 77 and the local diffusion coefficient D. By fitting the TeV 
data one can thus attempt to constrain, within the uncertainties given by the errors on 
the other measured quantities (namely, E$n, tage, M c i, and d) and by the assumptions 
made (e.g. the CR injection spectrum is assumed to be E~ 2 , while the energy dependence 
of D is assumed to scale as a power law of index 6 — 0.5), a combination of these two 
parameters (namely 77/D 3 / 2 ). The fact that the MCs have to be located within a distance 
Rd from the SNR can be verified a posteriori, and their exact location (unknown due to 
projection effects) can be tuned to match also the observed GeV emission. Given all the 
uncertainties above, our results have to be interpreted as a proof of concept of the fact 
that gamma ray observations of SNR/MC associations can serve as tools to estimate the 
CR diffusion coefficient. More detection of SNR/MC associations are needed in order to 
check whether the scenario described here applies to a whole class of objects and not 
only to a test-case as W28. Future observations from the Cherenkov Telescope Array will 
most likely solve this issue. 

Fig.[TJ from Gabici et al. 2010, shows a fit to the gamma-ray data for the three massive 
MCs in the W28 region. A simultaneous fit to all the three MCs is obtained by fixing 
a value for rj/D 3 / 2 , which implies that the diffusion coefficient of particle with energy 
3 TeV (these are the particles that produce most of the emission observed by HESS) is: 

2/3 



D(3 TeV) « 5 x 10 27 (£-\ cm 2 /s . (2.1) 



This value is significantly smaller (more than an order of magnitude) than the one nor- 
mally adopted to describe the diffusion of ~ TeV CRs in the galactic disk, which is 
w 10 29 cm 2 /s. This result remains valid (i.e. a suppression of the diffusion coefficient is 
indeed needed to fit data) even if a different value of the parameter 5 is assumed, within 
the range 0.3. ..0.7 compatible with CR data. 

The reason of this discrepancy between the average CR diffusion coefficient in the 
Galaxy and the one found in the vicinity of a SNR needs to be explained. A possible way 
to interpret these observations is given in the next Section. 
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Figure 1. Broad band fit to the gamma ray emission detected by FERMI and HESS from 
the sources HESS Jf80f-233, HESS Jf800-240 A and B (left to right), that coincide with three 
massive MCs. Dashed lines represent the contribution to the gamma-ray emission from CRs 
that escaped W28, dotted lines show the contribution from the CR sea, and solid lines the total 
emission. Distances to the SNR centre are 12, 65, and 32 pc (left to right). FERMI and HESS 
data points are plotted in black. No GeV emission has been detected from HESS J1800-240 A. 



3. Anisotropic diffusion of runaway cosmic rays 

Most of the studies aimed at predicting the gamma-ray emission from runaway CRs 
rely on the assumption of isotropic diffusion (see references in Sec. 2 and, e.g. Lee et al. 
2008, Torres et al. 2008, Rodrguez Marrero et al. 2008, Torres et al. 2010, Ellison & 
Bykov 2011, Li & Chen 2012, Ellison et al. 22012, Telezhinsky et al. 2012). However, 
the validity of the assumption of isotropic diffusion of CRs, adopted in the previous 
section, needs to be discussed. In fact, if the intensity of the turbulent field SB on scales 
resonant with the Larmor radius of particles is significantly smaller than the mean large 
scale field Bq (i.e. if 5B/Bq <C 1), then cosmic ray diffusion becomes anisotropic, with 
particles diffusing preferentially along the magnetic field lines (e.g. Casse et al. 2002 and 
references therein) . In the limiting (but still reasonable) case in which the perpendicular 
diffusion coefficient can be set equal to zero, the transport of CRs across the mean field 
is mainly due to the wandering of magnetic field lines (Jokipii & Parker 1969). 

To give a qualitative idea of the role that anisotropic diffusion can play in the studies 
of the CRs that escaped SNRs, let us consider an idealized case in which the escaping 
particles diffuse along a magnetic flux tube characterized by a very long coherence length 
(i.e. the magnetic flux tube is preserved for a long distance). In this case, after a time 
t particle will diffuse up to a distance Rd ~ y/2 D\\ x t along the tube (here Dm is the 
parallel diffusion coefficient of cosmic rays, not to be confused with the isotropic diffu- 
sion coefficient D adopted in the previous Section), while their transverse distribution 
will be equal to the radius of the SNR shock at the time of their escape, R s h, which is 
of the order of f=a 1-10 pc. Thus, the CR density in the flux tube will be proportional 
to ucr oc {RdR^h)" 1 instead of oc R^ 3 as in the isotropic case (see previous Section). 
It is easy to see that the estimates of the diffusion coefficient based on the two op- 
posite assumptions of isotropic and one-dimensional diffusion will differ by a factor of 
~ (Rd/Rsh) 4 ^ 3 , which can be much larger than an order of magnitude! Thus, it is of 
paramount importance to investigate how the interpretation of gamma-ray observations 
depends on the assumptions made concerning CR diffusion. 

As a first step, let us compare in Fig. [2] the results that are obtained if an isotropic 
diffusion coefficient is assumed, with the ones obtained for the anisotropic diffusion model 
considered in Nava & Gabici (2013). In both panels of Fig. [21 the SNR is located at the 
centre of the field and the color code refers to the excess of CRs with respect to the CR 
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Figure 2. Cosmic ray over-density above the galactic background around a supernova remnant 
(located at the centre of the panels) . A particle energy of E — 1 TeV and a time t = 10 kyr after 
the explosion are considered. The left panel refers to an isotropic diffusion coefficient of cosmic 
rays equal to D = 5 x 10 26 (-E/10 GeV) ' 5 cm 2 /s, while the right panel refers to an anisotropic 
diffusion scenario with Dm = 10 28 (-E/10 GeV) ' 5 cm 2 /s. The black cross marks the a posit ion at 
which the CR over-density is equal in the two panels. Figure from|Nava fc Gabici (2013)| 



sea. Over-densities are plotted for a particle energy of ITeV and for a time t = 10 kyr 
after the supernova explosion (see Nava & Gabici 2013 for more details on the model). 

The spatial distribution of CRs is strikingly different in the two scenarios: spherically 
symmetric in the left panel, and strongly elongated in the direction of the magnetic 
field flux tube in the right panel. A filamentary diffusion of CRs was also found in the 



numerical simulations by |Giacinti et al. (2012) The same parameters have been used to 
compute the over-densities in the two scenarios in Fig. [2j with the exception of the CR 
diffusion coefficient, which in the left panel has been assumed to be isotropic and equal to 
D(l TeV) w5x 10 27 cm 2 /s, while in the right one is assumed to be strictly parallel (i.e. 
CRs diffuse only along field lines) and equal to Dm w 10 29 cm 2 /s (a value similar to the 
average CR diffusion coefficient in the Galaxy). The choice of two significantly different 
values for the diffusion coefficients, with D <C D» has been made in order to obtain 
the same level of CR over-density in the vicinity of the SNR. As an example, the black 
cross in Fig. [2] identifies a position, located 30 pc away from the centre of the explosion, 
where the CR over-density is identical in the two panels. To get comparable values for 
the CR over-density, a much smaller (isotropic) diffusion coefficient D is needed in order 
to compensate for the larger solid angle over which CRs can propagate. This fact must 
be taken into account when interpreting the gamma-ray observations of molecular clouds 
illuminated by CRs escaping from SNRs. For example, a fit to the gamma-ray data from 
the MCs in the W28 regions has been obtained by Nava & Gabici (2013) by assuming a 
large diffusion coefficient of Dh (1 TeV) s=s 10 29 cm 2 /s, much larger than the (isotropic) 
one adopted in the previous section (see Fig. [I]). To conclude, the hint for a suppression of 
the diffusion coefficient in the vicinity of the SNR W28 obtained in the previous Section 
might depend on the assumption of isotropy of diffusion. If an anisotropic diffusion is 
adopted, a much larger diffusion coefficient can be assumed to fit gamma-ray data. 
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4. Conclusions and future perspectives 

I have shown how gamma-ray observations of MCs located close to SNRs can serve 
to support the idea that SNRs are the sources of CRs. Information on the CR diffusion 
coefficient can also be extracted from such observations, though the conclusions of these 
studies strongly depend on the (still unknown) isotropic or anisotropic nature of diffusion. 
To date, only two SNRs show gamma-ray emission clearly coming from outside the SNR 
shell: W28 (see above) and W44 (Uchiyama et al. 2012). So, further observations are 
needed in order to obtain solid constraints on the CR diffusion coefficient. Future facilities 
as the Cherenkov Telescope Array will play a key role in this direction. 

Theoretical studies are also needed in order to understand the details of CR propa- 
gation close to their sources and interpret correctly the gamma-ray observations. The 
diffusion of CRs along the magnetic field lines is most likely a nonlinear process, where 
the CRs themselves generate the magnetic turbulence needed to confine them. This effect 
is expected to be stronger in the vicinity of CR sources, due to the enhanced intensity 
of CRs. Pioneering works on the non-linear propagation of runaway CRs can be found 



in|Skilling (1970) and Hartquist & Morfill (1994) These studies have been revived by 



the recent results obtained from the gamma-ray observations of SNR/MC associations 
(see e.g. Ptuskin et al. 2008, Malkov et al. 2013), and promise to become one of the most 
important developments in this field. 
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